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ABSTRACT 

We have tested the appropriateness of two-sohton analytic metric to describe the 
exterior of all types of neutron stars, no matter what their equation of state or rotation 
rate is. The particular analytic solution of the vaccuum Einstein equations proved quite 
O . adjustable to mimic the metric functions of all numerically constructed neutron-star 

models that we used as a testbed. The neutron-star models covered a wide range of 
5^ I stiffness, with regard to the equation of state of their interior, and all rotation rates up 

{DJO to the maximum possible rotation rate allowed for each such star. Apart of the metric 

functions themselves, we have compared the radius of the innermost stable circular 
orbit i?iscO: the orbital frequency = ^ of circular geodesies, and their epicyclic 
frequencies fip, f2z, as well as the change of the energy of circular orbits per logarithmic 
change of orbital frequency AE. All these quantities, calculated by means of the two- 
soliton analytic metric, fitted with good accuracy the corresponding numerical ones 
as in previous analogous comparisons (although previous attempts were restricted 
to neutron star models with either high or low rotation rates). We believe that this 
particular analytic solution could be considered as an analytic faithful representation 
of the gravitation field of any rotating neutron star with such accuracy, that one could 
explore the interior structure of a neutron star by using this space-time to interpret 
observations of astrophysical processes that take place around it. 

Key words: gravitation - stars: neutron - equation of state - methods: analytical - 
relativity - accretion. 



1 INTRODUCTION 

The amount and accuracy of modern observations in var- 
ious parts of the electromagnetic spectrum has increased 
dramatically. In order to give astrophysically plausible ex- 
planations of the various problems related to the observa- 
tions we have to rely on theoretical assumptions that are at 
least as accurate as the data we are try ing to analyze . Ther e 
is a large class of observations (see e.g.. Ivan der Klid (|2006l ')') 
that are related to the astrophysical environment of com- 
pact relativistic objects (AGNs, LMXBs, etc). Furthermore, 
the anticipated successful gravitational wave detection will 
open a new window to observe such objects. In order to un- 
derstand these phenomena, one has to have a sufficiently 
accurate analytic description of the spacetime around such 
compact objects. If the central object is a black hole, there 
is a unique choice in the framework of general relativity: the 
Kerr spacetime. On the other hand the geometry around a 

* E-mail: georgios.pappasOguest.uni-tuebingen.de 
f E-mail: thapostolOphys. uoa.gr 



rotating neutron star is much more complicated, since it de- 
pends on many parameters related to the internal structure 
of the neutron star and the way it rotates. 

The assumption that the geometry around such an ob- 
ject is approximately th at of a Schwarzsch ild, or a Kerr 
metric (see for example Ivan der Khi l|2006D ') IS very sim- 
plistic and it may lead to erroneous conclusions about 
the actual astrophysical processes that ta ke place in th e 
close neighborhood of the star itse l f (cf. iPappasI (|2012f ) ; 
iPachon. Rueda fc Valenzuela-Toledol (|2012M . 

One can alternatively rely on numerical codes that are 
able to describe the geometry around a realistic neutron star 
in a tabular form on a given grid wi th sufficiently high accu- 
racy. T here are various groups (see IStergioulas fc FriedmanI 
(Il995|), and for an extended list of numerical schemes see 
Ste rgioulaj l|2003l ')'). that have acquired expertise in building 
relativistic models of astrophysical objects with adjustable 
physical characteristics and constructing the metric inside 
and outside such objects by solving numerically the full Ein- 
stein equations in stationary, ajdsymmetric cases. 
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Although studying astrophysical phenomena in a geo- 
metric background that has been constructed numerically 
is plausible, there are certain drawbacks in using such met- 
rics: (i) Computing various physical quantities of a system, 
like the orbital frequencies, or the innermost circular orbit, 
from a metric that is given in a tabular form is not very 
practical and is often plagued by numerical errors, (ii) As- 
trophysical observations from the environment of a compact 
object could be used to read the physical parameters that 
are related to the structure of the compact object such as its 
mass, its equation of state, its rotation, or to obtain the law 
of its differential rotation, etc. This would be very difficult 
to achieve with a metric that is numerically constructed. 
Solving the inverse problem by means of a numerical metric 
is a blind process that can not be easily led by physical in- 
sight. Instead an analytic expression for the corresponding 
spacetime would be much more preferable. 

There are various analytic metrics that have been used 
in the past to describe the exterior geometry of a neutron 
star. As mentioned above, the Schwarzschild metric is not 
accurate enough for rotating neutron stars, while the Kerr 
metric is good only for a collapsed object (a black hole) 
but it fails to describe the exterior of a neutron star, as 
comparisons of Kerr with numerical geome tries of rotating 
neutron stars bv lBerti fc StergioulasI (l2004h have shown. 

The Hartle-Thorne metric of iHartle fc Thornd l|l968h . 
which has been constructed as an approximate solution of 
the vacuum Einstein equations (VEE) for the exterior of a 
slowly rotating star, has been extensively used by various 
authors to de scribe neutr o n star s of low rotation rate (see 
for example iBerti et all (j2004l )). Finally, various other 
analytic solutions of VEE have been constructed and some 
of them have been used, especially during the last decade, 
to describe the ext er ior ge o metries of neutr on s ta rs (se e 
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Teichmiiller. Frob fc Maucheil ( 201ll)). Such solut ions are 



based on the formalism developed bv lErnsd (|l968al lbh which 
reformulates Einstein equations in the case of axisymmetric, 
stationary space-times . Manko et al. and SibgatuUin ( see th e 
articles of ISibeatulinl jll984h: iManko fc Sibeatuiinl lll993l); 
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I I995I ): iManko. Mielke fc Sanabria-Gomed ()2000l )') ha^ 

used various analytic methods to produce such space-times 
parameterized by various parameters that have different 
physical context depending on the type of each solution. 

Such an analytic solution, with its parameters appropri- 
ately adjusted to match numerical models of neutron stars, 
could then be used to describe the stationary properties of 
the space-time around the neutron star itself; that is, study 
the geodesies in the exterior of the neutron star. More specif- 
ically, from the analytic solution we could obtain bounds 
of motion for test particles orbiting the neutron star, find 
the location of the innermost stable circular orbit (ISCO), 
compute the orbital frequency of the circular orbits on the 
equatorial plane as well as the epicyclic frequencies around it 
and perform an y sort of dynamical analysis on the geodesies 
(for example see lLukes^G crakop oulos (2012 )). These proper- 
ties of the space-time could be used to study quantitatively 
astrophysical phenomena that take place in the vicinity of 
neutron stars, such as accretion discs. Inversely one could use 
the astrophysical observations related to such phenomena to 



determine the parameters describing the analytic space-time 
and from that acquire information for the central object. 

The central issue with analytic metrics is whether one 
can find solutions that are able to describe with sufficient 
faithfulness all kinds of rotating neutron stars; either slowly 
or rapidly rotating ones, or even differentially rotating ones. 

One solution that has been recently used by 
IStute fc Camenzindl (I2OO2I ) and later bv lBerti fc Stergioulad 
to describe the exterior space-time of rotat- 
ing neutron stars is the three - param eter solution of 
iManko. Mielke fc Sanabria-Gomea (|2000l ) (also mentioned 
as Manko et al.). Although this solution was shown to match 
quite well the space-time of highly rotating neutron stars, 
it failed to match the slowly rotating ones. The reason for 
failing to describe slow rotation is that in the zero angular 
momentum limit, this particular solution has a non van- 
ishing quadrupole moment, while one would expect slowly 
rotating neutron stars to be approximately spherically sym- 
metric. This problem of the Manko et al. solution was not 
considered disappointing by Berti and Stergioulas, since the 
space-time around slowly rotating stars could be described 
equally well by the Hartle-Thorne approximation. 

The three-parameter solution of Manko et al. is a spe- 
cial case of the so called two-soliton s olutio n, which was con- 
structed bv lManko. Martin fc Ruiz! (|l995l ). The two-soliton 
is a four-parameter analytic metric which, contrary to the 
previous one, can be continuously reduced to a Scwarzschild 
or a Kerr metric while it does not suffer from the problem- 
atic constraints of the Manko et al. solution with respect to 
the anomalous behavior of its quadrupole moment. Actually 
the first four multiple moments of the two-soliton solution 
can be freely chosen. Of course the analytic form of the 
two-soliton solution is not as compact as the Manko et al. 
solution, but this is the price one has to pay in order to cover 
the whole range of the physical parameters of a neutron star 
with a single analytic metric. 

In the present work, that constitutes the exten- 
sion an d com pletion of preliminary results presented by 
IPappasI l|2009[ ). that were suitably corrected with respect 
to the right extraction of the multipole moments of the 
numerical space-time as w as recently demonstrated by 
IPappas fc Apostolatod l|2012l) . we are using this two-soliton 
solution to describe the space-time around a wide range 
of numerically constructed rotating neutron stars. We use 
the numerical multipole moments to set the multipole mo- 
ments of the analytic space-time. Then we examine how 
well the two metrics match each other. Moreover we have 
performed comparisons between astrophysicaly relevant ge- 
ometric quantities produced from the numerical and the an- 
alytic space-times, like the position of the innermost stable 
circular orbit (ISCO), the orbital frequencies, the epicyclic 
frequencies that are related to slightly non-circular and 
slightly non-equatorial orbits, and the change of energy of 
the circular orbits per logarithmic change of the orbital fre- 
quency, /S.E. The overall picture is that the new metric 
matches the numerical one with excellent accuracy for all 
rotation rates and all equations of state. 

The rest of the paper is organized as follows: In Sec. 2 
the proposed analytic solution (two-soliton) is briefiy pre- 
sented and some of its properties are thoroughly analyzed. 
The parameter space of the two-soliton is investigated and it 
is shown how to obtain the limiting cases of Schwarzschild, 
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Kerr and Manko et al. A brief discussion of the physical 
properties of the space-time such as the presence of singu- 
larities, horizons, ergoregions and regions of closed timelike 
curves (CTCs) is also given. In Sec. 3 we show how we match 
the analytic two-soliton solution to a specific numerical one 
by matching the first four multipole moments, and show why 
this is generally the best choice. In Sec. 4 we discuss various 
criteria that could be used to compare the two metrics. Fi- 
nally, in Sec. 5 the final comparison criteria and the results 
of the corresponding comparisons are presented. In Sec. 6 we 
give an overview of the conclusions obtained by our study. 



2 THE TWO-SOLITON SOLUTION 

The vacuum region of a stationary and axially symmetric 
space-time can be descr ibed by the Papape trou line element, 
which was first used bv lPapapetroul (|l953f ). 



-/ {dt - u}d4>f + \e-^' {dp^ dz^) 



(1) 



where /, tj, and 7 are functions of the Weyl-Papapetrou 
coordinates (p, z). By intr oducing the complex potential 
E{p, z) = /(p, z) + til^(p, z). lErnstI 1 1968a') reformulated the 
Einstein field equations for this type of space-times in a con- 
cise complex equation 



{Re{£))\/^£ = V£ ■ VS. 



(2) 



The real part of the Ernst potential £ is the metric 
function /, which is also the norm of the timelike Killing 
vector related with stationarity, while 1/) is a scalar potential 
related with the twist of the same vector. 

A general procedure for generati ng solutiori s of th e 



Ernst equations w a s de v eloped by SibgatulinI 
_Manko & Sibgatulin' ('l993''): jRuiz. Manko fc Martini 
Manko. Martin 



1984 ): 



199511: 

Ruiz, ((iggsl ). Each solution of the Ernst 



equation is produced from a choice of the Ernst potential 
along the axis of symmetry of the metric in the form of a 
rational function 



£{p = 0,z) = e{z) = 



Piz) 
R{zy 



(3) 



where P{z), R{z) are polynomials of z of order n with com- 
plex coefficients in general. The algorithm developed by 
I Ruiz. Manko fc Martini (Il995l ) works as follows: First the 
Ernst potential along the axis is expressed in the form 



(4) 



where /3fc are the roots of the polynomial R{z) and ek are 
complex coefficients appropriately chosen so that the lat- 
ter form of e{z) (eq. H} is equal to the former one (eq. |3]). 
Subsequently one determines the 2n roots of equation 



e{z) + e'iz) = 0, 



(5) 



where * denotes complex conjugation. These roots are de- 
noted as ^fe, with k — 1,2, . . . ,2n a nd from these o ne defines 
the 2n complex functions Rk = ^ p^ -\- (z — ^k)^- AH these 
functions and roots are then used as building blocks for the 
following determinants 



E± = 



G - 



1 1 

±1 

±1 






Rl+^1- Z ■■■ R2n+^2n~ Z 
-1 



(6) 



H = 





z 1 

-Pi 

el 



(7) 



(8) 



€l-/3r. 



C2„-/3l 



(9) 



where C is the 2n x 2n matrix 



Hi 



C = 



R2r, 



e2„-/3 



i2„ 



(10) 



e2„ 



The Ernst potential and the metric functions are finally ex- 
pressed in terms of the determinants given above as: 

E-\-E/Z- + Ej^E— 



£{p,z) 



uj{p,z) 



2E-E*_ 
E+E*_ + EXE- 
2KoK*Y{ll^Rk 

2Q[E-H* - EIG] 



(11) 
(12) 

(13) 

(14) 



E+El + E^E^ 

We should note that due to the form of the metric func- 
tions, the parameters ek and their complex conjugates el 
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that appear in the determinants cancel out (the O^^j Cfee^. 
is a common factor of all products of determinants that show 
up in the metric functions), so they do not affect the final 
expressions. 

The vacuum tw o -solito n solution (proposed by 
iManko. Martin fc Rui j (jl995l 'l) is a special case of the 
previous general axisym metric solution that is obtained 
from the ansatz (see also ISotiriou fc Pappad (|2005l )) 



{z- M - ia){z + ib)-k 



(15) 



[z + M - ia){z + ib)- k' 

where all the parameters M, a, k, b are real. From the Ernst 
potential along the axis one can compute the mass and mass- 
current moments of this space-time. Particularly, for the 
two-soliton space-time the first five non-vanishing moments 
are: 



Mo 


= M, 


Ma = 


-(a^ -A:)M, 




Mi 




4 

a - 


- {3a^ - 


- 2ab + b^)k + k^ - ifcM^ 
7 


M 


Ji 


= aM, 


Ji = 


-[a^ - {2a-b)k]M. 


(16) 



The mass moments of odd order and the mass-current mo- 
ments of even order are zero due to reflection symmetry 
with respect to the equatorial plane {z = 0) of the space- 
time (this is actually ensured by restricting all parameters 
of eq. (|15[l to assume real values). From the moments we see 
that the parameter M corresponds to the mass monopole 
of the space-time, the parameter a is the reduced angular 
momentum, k is the deviation of the reduced quadrupole 
from the corresponding Kerr quadrupole (the one that has 
the same M and a), and b is associated with the deviation 
of the current octupole moment from the current octupole 
of the corresponding Kerr. 

For the two-soliton ansatz (|15|) . the characteristic equa- 
tion ((5]) takes the form. 



(M^ 



6^ + 2k)z^ + (fc - abf 



0. (17) 



Since the coefficients of the polynomial are real, the roots 
can be either real, or conjugate pairs. The four roots of p7p 
can be written as. 



where, 



with 



K± = y^M^ ~a^-b'^ + 2k± 2d, 
and 



(18) 



(19) 



(20) 



d = y/{k- ab)^ -b'^M'^. (21) 

Using these symbols for the four roots we redefine the four 
corresponding functions Rk as 



R± = v/p2 + (z±C+)2, r± = ^p2 + {z±^.y. 



(22) 



Next we proceed to classify the various types of solu- 
tions depending on whether the four roots have real, purely 
imaginary or complex values. This classification is outlined 
in Table 1. 

• Case I: This case is characterized by two real roots 



The Kerr family of solutions, which corresponds to fc = is 
definitely not included in this family of solutions. 

la: This subfamily of case I is the simplest to compute, 
since all functions R± , r± are real. 

lb: This is a degenerate case where the roots ^± coin- 
cide. The degeneracy is due to k_ being zero which corre- 
sponds to the parameter constraint M^— a^— 2fc— 2d = 
0. In such degenerate cases, the computation of the met- 
ric function is not straightforward since the expressions 
for the metric become indeterminate, and a limiting pro- 
cedure should then be applied. In the reduced-parameter 
space {a/M,b/M, k/ M^), the previous constraint corre- 
sponds to a two-dimensional surface. 

• Case II: In this case, the roots ^± are either complex or 
imaginary, since k'L < 0. Furthermore this means that there 
are nonvanishing values of (p, z) where the functions R±,r± 
assume zero value, which then leads to singularities at the 
corresponding points. 

Ila: This sub-case, as with Case I, belongs to a class of 
solutions that cannot have a vanishing parameter fc. 

lib: Here, a degeneracy shows up again as in Case lb, 
which admits the same treatment (limiting procedure) as 
in the former situation. Contrary to all previous cases. 
Case lib admits a Kerr solution that belongs to the hyper- 
extreme branch (|a| > A/). Similarly to Case lb this so- 
lution is also represented by a two dimensional surface in 
the reduced-parameter space. 

lie: This case is similar to the previous one, without the 
degeneracy in the roots ^± . It also includes hyper-extreme 
Kerr solutions. 

• Case III: In this case one of the C+iC- is real while the 
other one is imaginary. Thus the same type of singularity 
issues, as in Case II arise. In particular such problematic 
behavior shows up on the z — plane. The Kerr and the 
Schwarzschild solutions lie entirely within this family of so- 
lutions. 

• Case IV: All types of solutions belonging to this case 
are degenerate (there is a special constraint between the 
parameters) and as such are probably of no interest to re- 
alistic neutron stars. Subcases IVa and IVc have one dou- 
ble root (^_ — 0), while subcase IVb has a quadruple root 
(^-i_ = ^_ = 0) and the computation of the metric functions 
needs special treatment. We should also note that cases IVb 
and IVc include the extreme Kerr solution (|a| = M) as a 
special case. 

As we can see from the classification, the two-soliton 
solution can produce a very rich family of analytic solutions 
with the classical solutions of Schwarzschild and Kerr be- 
ing special cases of the general solution. Also the Manko 
et al. solution of lManko. Mielkc & Sanabria-Gomcz ( 200^ ) 
that has been used p r eviou sly bv lBerti fc Stergioulad (|20o4 ): 
IStute fc Camenzindl (|2002l l to match the exterior space-time 
of rotating neutron stars is a special case of the two-soliton 
solution as we will see next. 

All types of solutions discussed above can be repre- 
sented in a three dimensional parameter space, the reduced- 
parameter space that was mentioned in Case lib. Although 
the two-soliton solution is characterized by four parame- 
ters, one of them, the monopole mass M, is simply a seal- 
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Table 1. Classification of the various two-soliton solutions depending on the values of the 
parameters d, k±, R± and r±. The table also shows the various conjugation relations 
between the parameters. The types of solutions that have degeneracies are indicated with an 
asterisk (*). 3? means real, 9? means imaginary, and C means complex. 
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Figure 1. The first two plots show the parameter space (a,b,k) of the two soliton solution for a particular mass (M = 1) from two 
different view points so that the intricate foldings of the surface are better understood. The 2D surface plotted is the constraint of k 
which corresponds to the solution of Manko et al. used by Berti & Stergioulas (2004). The third plot is the k = plane of the parameter 
space which corresponds to all Kerr solutions. The plots clearly show that the Manko et al. solution has no set of parameters to describe 
the case k = which corresponds to the Kerr and the Schwarzschild (for a = 0) solutions, since there is no intersection of the constraint 
of k and the Kerr plane in the appropriate range of parameters. The two hyperbolas plotted on the plane k = axe the only points were 
the constraint of k touches the plane tangentially. As we can see these hyperbolas correspond to \a\ > M, i.e., to hyper-extreme Kerr 
spa<;e-times. 
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ing parameter which can be used to reduce the rest of the 
parameters to dimensionless ones. The three dimensionless 
parameters thus formed, {a/M,h/M,k/ M"^), are related to 



the multipole moments (see eq. (|16|l 'l of the corresponding 
space-time in the foUowing way: The first parameter a/M 
is the spin parameter (where a is the reduced angular mo- 
mentum) which is the only parameter, besides the mass, 
that uniquely characterizes a Kerr space-time. The second 
parameter k/M"^ expresses the deviation of the quadrupole 
moment of the solution from the quadrupole moment of the 
corresponding Kerr (the one with the same a/M value); an 
increase of the value of kjM^ produces solutions that are 
less oblate than Kerr. The final parameter hjM controls in 
a linear fashion the current octupole moment. The actual de- 
viation of the two soliton octupole moment from the Kerr oc- 
tupole moment depends on all three parameters a/M, k/M^ 
and b/M. Of course the higher moments are also affected by 
these parameters. 

In this three dimensional parameter space, the plane 
k/M"^ — corresponds to all types of Kerr solutions. This 
is clear from the form of the Ernst potential along the axis, 
where if one sets = it reduces to the Ernst potential of 
the Kerr solution, 



e{z) 



M ■ 



(23) 



z + M - ia' 

Obviously in this case the parameter b/A'I is redundant; thus 
each line a/M = const, which is parallel to the b/AI axis on 
the plane k/M — 0, corresponds to a single Kerr (modulo 
the mass of the black hole) . 

As mentioned in the introduction , the solution of 
iManko. Mielke fc Sanabria-GonieS (|2000l ') has been used to 
describe the exterior of rotating neutron stars. As it was 
briefly discussed above this solution is included in the two- 
soliton solution and can be obtained by imposing a specific 
constraint on the two-soliton parameters. The Manko et al. 
solution is obtained by setting 

^-bf]-jj^^:-^+ab. (24) 



M'' 



(a-fe)2 

This constraint defines a surface in the three parameter 
space {a/M, b/M, k/M^) (see Figure[T]). The particular solu- 
tion, depending on the values of a, b, falls under either Case 
lb or Case lib, where either k_ or k+ is equal to zero, re- 
spectively. We should note that the Manko et al. solution 
is the union of these two cases. By substituting the above 
expression for k (eq. ((24])) in the formula for the quadrupole 
moment (I16p . the quadrupole moment takes the following 
value when a = 0, 



Ma 



M (M^ + b^ 



(25) 



4 Af2 - b2 ■ 

This is why the quadrupole moment of the Manko et al. so- 
lution does not vanish in the limit of zero rotation. From the 
above expression one can see that the metric is not spher- 
ically symmetric as one would expect for a non rotating 
object. Especially for |6| < M the metric is oblate while for 
|6| > M the metric is prolate. 

This anomalous behavior of the quadrupole moment is 
an important drawback for using the Manko et al. solution to 
de scribe every rotating neutr on star and it was pointed out 
bv lBerti fc StergioulasI (|2004l V In fact this analytic metric is 
good to match only rapidly rotating neutron stars. 



As shown bv lManko. Mielke fc Sanabria-Gomed (|2000l ) 
this particular metric turns into a Kerr metric if b^ = a^ — 
M^. Since all the parameters are assumed real, then this 
corresponds to a hyperextreme Kerr metric (a ^ M). In 
Figure [T] this is represented by the two hyperbolas that lay 
outside the strip \a\ < M on the Kerr plane (the plane k = 0) 
along which the intricate surface of the Manko et al. solution 
tangentially touches the corresponding plane. 

The two-soliton solution, which we will thoroughly 
study later on, is a much better metric to describe the ex- 
terior of an arbitrary rotating neutron star than the Manko 
et al. solution because (i) the former has 4 independent pa- 
rameters (compared to the 3 independent parameters of the 
latter one) that offer more flexibility to adjust the metric, 
and (ii) these 4 parameters are able to cover the whole space 
of the flrst 4 moments of the space-time, while the first 4 mo- 
ments of the latter metric are actually correlated with each 
other through the dependence of k on the 3 independent 
parameters M, a, b, that was mentioned previously. Exactly 
this restriction renders the Manko et al. solution inappro- 
priate to describe slowly rotating neutron stars. 

Before closing this section we will give a brief descrip- 
tion of the space-time characteristics of the two-soliton so- 
lution for the range of parameter that we are going to use. 

A horizon of a space-time is the boundary between the 
region where stationary observers can exist and the region 
where such observers cannot exist. For a stationary and 
axially symmetric space-time, the stationary observers are 
those that have a four-velocity that is a linear combination 
of the timelike and the spacelike Killing vectors that the 
space-time possesses, i.e. 



(26) 



where ^'^ , rj^ are the timelike and spacelike Killing fields 
respectively and is the observer's angular velocity. The 
factor 7 is meant to normalize the four-velocity so that 
g^vU^u" ~ —1. In order for the four- velocity to be time- 
like, 7 should satisfy the equation 



(27) 



and it should be 7 > 0, which corresponds to an Q, taking 
values between the two roots 



0,4 



gttc 



g<t><i> 



(28) 



This condition can not be satisfied when (gttp)^ — 
5tt50</) ^ 0. Thus the condition (gt^)^ — gttg44> = defines 
the horizon. In the case of the two-soliton, expressed in the 
Weyl-Papapetrou coordinates, this condition corresponds to 
having p — 0, since = {gtct/f — gttg<i>4,- Thus the issue of 
horizons is something that we will not have to face; in these 
coordinates the whole space described corresponds to the 
exterior of any possible horizon. 

Another issue is the existence of singularities. Singular- 
ities might arise where the metric functions have infinities. 
From equations (|ll)ll4p one can see that singularities might 
exist where the functions R± = p^ + {z ± , r± — 
p^ + {z ± C-)^ go to zero, or where the determinant E- 
goes to zero, or where E+E*_ -\-E'^E- goes to zero. Whether 
or not some of these quantities vanish depends on which 
Case the solution belongs to. A thorough investigation of 
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Figure 2. Typical types of tlie surfaces that define tlie static limit and the regions of CTCs for the various cases of the two-soliton. On 
the left panels we have two typical figures that correspond to Case la, in the middle panels the figures correspond to Case Ila and on 
the right panels the figures correspond to Case III. The solid curves correspond to the static limit and the dotted curves correspond to 
the boundary of the regions with CTCs. For the Case la, the region of CTCs is in contact with the axis of symmetry while the same 
also applies to the surface of the static limit. For the Case Ila the static limit is in contact with the axis, but the region of CTCs has 
been detached. Finally for the Case III the inverse is true, that is the region of CTCs is in contact with the axis of symmetry while the 
static limit is detached from the axis. In all cases the upper plots correspond to slow rotation while the lower plots correspond to faster 
rotation. As the rotation rate increases the upper configurations are continuously deformed to the lower ones. 



the singularities of the two-soliton is out of the scope of this 
analysis; so we should only point out that for all the neutron 
star models that we have studied and the corresponding pa- 
rameters of the two-soliton solution, any such singularities, 
when present, are always confined in the region covered by 
the interior of the neutron star and thus they do not pose 
any computational problems in our analysis. 

The final issue is with regard to the existence of er- 
goregions (i.e., regions where g^v^'^^" ~ gtt < 0) and re- 
gions with closed timelike curves (CTCs) (i.e., regions where 
gabrj^'ri^ = g<i><t, < 0). In Figure[2]we have plotted the bound- 
ary surfaces of such regions for the two-soliton metric. One 
can see that there are three distinct topologies for these sur- 
faces that are observed for the different two-soliton Cases. 
In any case though, for all the models used here, these sur- 
faces are again confined at regions where the interior of the 
neutron star lays. 

In the following sections we analyze the method that we 
are going to use to obtain the right values for the parame- 
ters of the two-soliton solution for each neutron star model 
and compare its properties with the corresponding numeri- 
cal metric. In order to do that, we have constructed several 



sequences of numerical neu tron star models with th e aid o f 
the RNS numerical code of IStergioulas &: FriedmanI l|l995h . 
The numer ical neut ron star models used are the same mod- 
els used bv lPappas fc Apostolatos (2012|) for demonstrating 
how to correct the numerical multipole moments. They are 
produced using three equations of state (EOSs), i.e., AU, 
FPS and L. The scope of using these models is twofold. First 
we are using them to provide the appropriate parameters 
describing realistic neutron stars in order to build the cor- 
responding analytic metrics. Then we use them as a testbed 
against which to compare the analytic metrics and thus test 
their accuracy. As we have already mentioned, we will use as 
matching conditions between the analytic and the numeri- 
cal metrics the first four non-zero multipole moments. For 
the neutron star models that we have studied, the corre- 
sponding analytic space-times that are produced belong to 
three of the Cases of the aforementioned classification, i.e., 
to Cases la, Ila and III. 
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Figure 3. Contour plots that point out what is the best choice for the parameters of an analytic metric so that it matches well a 
numerical one. The left plot shows the contour plots of the overall mismatch aij between the analytic and the numerical metric for the 
tt (black curves) and t(f> components (grey curves), respectively, as a function of the fractional deviation of the quadrupole, 5M2, and 
the current octupole, SSg, of the analytic metric from those calculated directly from the numerical metric (assuming the same mass and 
angular momentum though). Since the contours of the att (almost horizontal) are orthogonal to the ones of (j^^ (almost vertical), the 
combination indicates an optimum choice for the multipole moments of the analytic space-time. That choice is the moments that have 
zero deviation from the moments of the numerical space-time. The next two contour plots are similar contour plots but the contours 
correspond to different quantities. The middle one is for the relative difference between the analytic and numerical -RiscOi while the right 
one is for the overall difference between the analytic and numerical orbital frequency fl (defined analogously to the overall mismatch of 
the metric components). Both these latter plots are consistent with the first one. All plots correspond to the model #15 for the AU EOS 
of the models presented by Pappas & Apostolatos (2012) and they give a representative picture of what's happening with all numerical 
models that we have computed. 



3 MATCHING THE ANALYTIC TO THE 
NUMERICAL SOLUTION 

When one attempts to match an analytic solution to a nu- 
merical one, it is desirable to find a suitable matching crite- 
rion that would be characteristic of the whole structure of 
the particular numerical space-time, instead of just a finite 
regio n of it. That i s, the matching should be global and not 
local. [Berti fc Stergioulas (2004) have argued that a suitable 
global condition should be the matching of the first few mul- 
tipole moments. Indeed, the full set of mu l tipole moment s 
(as defined relativistically b y lGerochl l|l970l ) ; iHanse ^ (119741) : 
iFodor. Honselaers fc Peried (|l989l )) of a stationary and axi- 
ally symmetric space-time can fully specify the Ernst poten- 
tial on the axis of symmetry. On the other hand, when the 
Ernst potential along the axis of symmetry is given, there is a 
space-time which is unambigu ously specified by tha t Erns t 
potential as it was sh own by IXanthopoulosI (|l979l . Il98ll ) ; 
iHauser &: ErnstI (|l98ll ). Thus, the full set of multipole mo- 
ments are uniquely characterizing a space-time and they can 
be used as a global matching condition. 

When the space-time of a neutron star model is 
constructed from a numerical algorithm, one can eval- 
uate its mass moments M,Q,... and current moments 
JjSs,... with an accuracy depending on the grid used 
to present the numeri c al m etric (for further discussion 
see [ Berti fc Stergioulas! (|2004l ). and iPappas fc Apostolato"^ 
Practically, the first few numerically evaluated mo- 
ments can be used as matching conditions to the analytic 
space-time. The first four nonzero multipole moments of the 
two-soliton solution as a function of its parameters M, a, b, k 
are shown in eq. (|16p from which it is clear that once we spec- 
ify the mass and the angular momentum of the space-time, 
the parameter k is uniquely determined by the quadrupole 



moment Q = M2 ,while the parameter h is uniquely de- 
termined by the current octupole Ss = J3. Thus, having 
constrained the four parameters of the two-soliton, we have 
completely specified an analytic space-time that could be 
used to describe the exterior of the particular neutron star 
model. What remains to be seen is how well do the prop- 
erties of the analytic space-time compare to those of the 
numerical one. 

At this point there is an issue that should be addressed. 
Having specified the first four non-zero moments of the two- 
soliton metric, we have fixed all the higher moments of the 
space-time in a specific manner related to the particular 
choice of the analytic metric. These higher moments will 
probably deviate from the ones of the numerical space-time. 
So the question is, could one make a better choice when 
trying to match the analytic to the numerical space-time 
than the one of setting the first four analytic moments ex- 
actly equal to the first four numerical moments? To an- 
swer this question we have performed the following test. 
For several numerical models of uniformly rotating neutron 
stars that we have constructed, we formed a set of two- 
soliton space-times for each neutron star model that have 
the same mass M and angular momentum J with the nu- 
merical model, while the quadrupoles and the current oc- 
tupoles of each single two-soliton space-time take the values 
M^"' = M2'"'(l - SM2) and 5^"' = S^l\l - SS^) respec- 
tively with various 5M2 and 53^ values. The quantities 5M2 
and 5S->, denote the fractional differences of the correspond- 
ing analytic moments of each two-soliton space-time from 
the numerical one. Then for each one of these sets of mo- 
ments we calculated the overall mismatch between the ana- 
lytic aB^_th«jiurn«^ which are defined 
(see lPappas fc Apostolato3 l|2012l )) as 
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{dij - 9ij) dr 



iRs 



(29) 



where Rs is the radius r at the surface of the star, and have 
thus constructed contour plots of atj on the plane of SS3 
and 6M2 , like the ones shown in Figure [31 The same type 
of contour plots were drawn for other quantities as well, like 
the relative difference of the i?isco and the overall difference 
between the analytic and numerical orbital frequency Q for 
circular, equatorial orbits (defined in the same fashion as it 
was defined the overall metric mismatch in eq. (|29p ). The 
contour plots for the particular neutron star model that is 
shown in Figure [3] are typical of the behavior that we ob- 
served in all models. An important result is that the contours 
for the overall mismatch att combined with the contours for 
at<i> give us the best choice for matching a numerical to an 
analytic space-time; namely, the equation of the first four 
multipole moments between the two space-times. That is 
because gtt seems to be sensitive mainly in deviations from 
the numerical quadrupole and thus the contours appear to 
be approximately horizontal and parallel to the axis of SS3, 
while gttj, seems to be sensitive mainly in deviations from the 
numerical current octupole and thus the contours appear to 
be approximately vertical and parallel to the axis of 5AI2 
(the contours of (Jt^ are almost orthogonal to the contours 
of att)- 

We should note here that the exact position of the op- 
timum point in the contour plots of att, atrj, did not deviate 
from (0, 0) by more than 3-4 per cent in all cases studied. 
The largest deviations showed up in some of the fastest ro- 
tating models. 

The conclusion is that what we expected to be true 
based on theoretical considerations turns out to be exactly 
the case after the implementation of the aforementioned test. 
Therefore in what follows, we will set the first four non-zero 
multipole moments of the analytic space-time equal to that 
of the corresponding numerical space-times. 



4 CRITERIA FOR THE COMPARISON OF 
THE ANALYTIC TO THE NUMERICAL 
SPACE-TIME 

Once we have constructed the analytic metric, appropriately 
matched to the corresponding numerical one, we proceed to 
thoroughly compare the two space-times. In order to do that, 
we should try again to use criteria that are characteristic 
of the geometric structure of the whole space-time, and if 
possible coordinate independent. It would be preferable if 
these criteria are also related to quantities that are relevant 
to astrophysical observations. Thus, if two space-times are in 
good agreement, with respect to these criteria, they could 
be considered more or less equivalent. On the other hand, 
such criteria, as well as possible observations associated to 
them, could be used to distinguish different space-times and 
consequently different compact objects that are the sources 
of these space-times. 

As a first criterion of comparing the metrics we will 
use the direct comparison between the analytic and the nu- 
merical metric components themselves. Although the metric 
components are quantities that are not coordinate indepen- 
dent, they have specific physical meanings and can be related 



to observable quantities. Thus, the gtt component is related 
to the gravitational redshift of a photon and the injection en- 
ergy of a particle. The gtif, component is related to the frame 
dragging effect, and the angular velocity uu = —gt<t>/ g<t,<t> of 
the zero angular momentum observers (ZAMOs). Finally the 
component is related to the circumference of a circle at 
a particular radial distance and defines the circumferential 
radius -Rcirc = C/27r = y/g^xj,- Also g^^ together with gr,- are 
used to measure surface areas. So, if the relative difference 
between the numerical and the analytic metric components 



gtt = = /oj. 



(30) 



is small, then one could consider the analytic metric as a 
good approximation of the numerical metric. 

Another criterion for comparing an analytic to a im- 
merical space-time is the location of the innermost stable 
circular orbit (ISCO). Particles moving on the equatorial 
pla ne are govern ed by the equation of motion (see for exam- 
ple |l99i) ) 



9pp 



= 1 - 



E^g^4, + 2ELgt4, L^gtt 



= V{p), (31) 



where E and L are the conserved energy and angular mo- 
mentum parallel to the axis of symmetry, per unit mass. 
V{p) is an effective potential for the radial motion and in 
the case of orbits that are circular, we additionally have the 
conditions dp/dr — and d^ p/dr^ = 0, which are equivalent 
to the conditions for a local extremum of the potential, i.e., 
V{p) = 0, and dV{p)/dp = 0. The radius of the ISCO is eval- 
uated if we further demand the constraint d^V{p)/dp^ — 0, 
the physical meaning of which is that the position of the 
circular orbit is also a turning point of the potential. From 
these three conditions we can evaluate a specific pisco and 
from that the -Risco = \/300(pisco), which we then com- 
pare to the corresponding numerical one. The position of the 
ISCO is of obvious astrophysical interest since it is the inner 
radius of an accretion disc and recently it has been used to 
evaluate the rotation parameter of black holes from fitting 
the continuo us spectrum of the accretion disc around them 
(see work bv lShafee etahl (|2006l )). 

Another criterion for comparing the metrics can be the 
orbital frequency of circular equatorial orbits Q. The orbital 
frequency is given by the equation 



n(p) 



'3*0, p + \/ {gt4,,py - gtt,pg4>4,,p 

9<t><t>.p 



(32) 



Apart from the orbital frequency one could also use the pre- 
cession frequencies of the almost circular and almost equa- 
torial orbits, i.e., the precession of the periastron Q,p and 
the precession of the orbital plane Q.^ . These frequencies are 
derived from the perturbation of the equation of motion 



9pp 



(33) 



around the circular equatorial orbits. In this expression, 
V{p, z) is the same effective potential which was defined 
in the second equation of (|3ip the z dependence of which 
now has not been omitted as in eq. (|3ip . The perturbation 
frequencies derived from the above equation are then given 
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with respect to the metric functions as 



5 RESULTS OF THE COMPARISON 



+ 



igt4> + g4>4'^f ( ^ 



gt<t, 



(34) 



where the index a takes either the value p or 2: to obtain 
the frequency of the radial or the vertical perturbation re- 
spectively. These expressions are evaluated on the equato- 
rial plane (at z = Q); thus they are functions of p alone. The 
corresponding precession frequencies are given by the dif- 
ference between the orbital frequency and the perturbation 
frequency: 



(35) 



These quantities are quite interesting with respect to astro- 
physical phenomena as well. More specifically, they can be 
associated to the orbital motion of material accreting onto 
a compact object through an accretion disc. These very fre- 
quencies have been proposed to be connected to the observed 
quasi-periodic modulation (QPOs) of the X-ray flux of ac- 
cretio n discs that are pres e nt in X-ray binar i es (se e 
' 2000[) : lvan der Klid (|2006l ): lBoutloukos et all (|2006l '): 
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Finally, the last criterion that we will use to compare 
metrics is the quantity AE of circular orbits, which expresses 
the energy difference of the orbits per logarithmic orbital 
frequency interval as one moves from one circular orbit to 
the next towards the central object. This quantity is defined 



(36) 



where the energy per unit mass E is given by the expression 

-gtt - gt4>^ 



E = 



6^2 



(37) 



The quantity AE is a measure of the energy that a parti- 
cle has to lose in order to move from one circular orbit to 
another closer to the central object so that the frequency in- 
creases by one e-fold. The quantity, A_E, is associated to the 
emissi on of gravitational radiation and was used by iRvanI 
||1995D to measure the multipole moments of the space-time 
from gravitational waves emitted by test particles orbiting 
in that background. The same quantity can also be associ- 
ated to accretion discs and in particular, in the case of thin 
discs, it would correspond to the amount of energy that the 
disc will radiate as a function of the radius from the central 
object and thus it will be related to the temperature profile 
of the disc and consequently to the tot al lumi n osity of the 
disc (for a review on accretion discs see iKroUkI jimi)'). 

The last set of criteria, i.e., the frequencies and AiJ, 
are related to specific observable properties of astrophysical 
systems, in particular of accretion discs around compact ob- 
jects; thus they are very useful and relevant to ast rophysics 
(for an application see the work bv lPappasI l|2012h l. 



The results of the comparison between the analytic and the 
numerical metrics, describing the exterior of realistic neu- 
tron stars, that are presented here is indicative of all com- 
parisons performed for every numerically constructed neu- 
tron star model. The models that we have used as a testbed 
of comparison are briefly presented in the Appendix an d are 
the same models used in lPappas fc Apostolatoj (|2012l ). 

For illustrative reasons we have also plotted compar- 
isons between the numerical and the Manko et al. solu- 
tion (iManko. Mielke fc Sanabria-G6rne3 l|2000l )) which was 
used by Berti fc StergioulasI 1 2004 ). as well as compar- 
isons between the numer ical and the Hartle-Thorne metric 
l|Hartle fc Thornd l|l968l )'). The reason for using these two 
metrics is that on the one hand the Hartle-Thorne metric 
is considered to be a good approximation of slowly rotat- 
ing relativistic stars and on the other h and the Manko et al. 
metric has been shown bv lBerti fc Stergioulas (2004) to be a 
good approximation for relativistic stars with fast rotation. 
In cases with slow rotation rates, for which corresponding 
models of the Manko et al. metric cannot be constructed, 
we have used only the Hartle-Thorne metric to compare, 
though. In the case of models with fast rotation besides the 
Manko et al. metric, we have used the Hartle-Thorne metric 
as well, for these cases we treated the Hartle-Thorne met- 
ric as a three parameter exterior metric, where the three 
parameters are the mass M , the angular momentum J and 
the reduced quadrupole q = M2/M^ of the neutron star. It 
should be noted that this is not a consistent way to use the 
Hartle-Thorne metric, since the quadrupole and the angular 
momentum in the Hartle-Thorne cannot take arbitrary val- 
ues, while the metric is essentially a two parameter solution 
(parameterized by the central density of the corresponding 
slowly rotating star and a small parameter e that corre- 
sponds to the fraction of the angular velocity of the star 
relative to the Keplerian angular velocity of the surface of 
the star) that has to be properly matched to an interior solu - 
tion following the procedure described bv lBerti et al.l |20o3). 
Here we are taking some leeway in using Hartle-Thorne for 
fast rotation since it is used simply for illustrative purposes 
and not in order to draw any conclusions from it. 

In Figure |4] we present the comparison of the various 
analytic metric functions (using the two-soliton, the Manko 
et al., and the Hartle-Thorne solutions) to the correspond- 
ing numerical ones for a single model constructed using the 
EOS AU (model #10 of the AU EOS the characteristics 
of which are presente d in T able II of the supplement of 
IPappas fc Apostolatos! l|2012l )). The figures display the rel- 
ative difference between the various analytic and the nu- 
merical metric functions gtt,gt<p, ^/gH> s-^d — 5pp the 
equatorial plane, as well as the function gu on the axis of 
symmetry. 

The general picture we get from these figures is typi- 
cal for all models constructed using all three equations of 
state, that is AU, FPS and L. The overall comparison of 
the two-soliton to the numerical metrics shows that this 
analytic metric is an excellent substitute of the numerical 
space-time both for slow and fast rotating models, with an 
accuracy that is everywhere outside the neutron star always 
better than about 1/1000 for all the metric functions (there 
is an exception to that for the comparison of the gu met- 



An analytic spacetime exterior to neutron stars 11 




1.0 1.5 2.0 2.5 3.0 3.5 1.0 1.5 2.0 2.5 3.0 

log(i-) (km) log(z) (km) 



Figure 4. These are the plots of the logarithm of the relative difference of the analytic to the numerical metric components gtt, gt4>, 
y/g^ = i?circi S'lid 9zz on the equatorial plane, as well as gtt on the axis of symmetry. The plots shown here are indicative of all 
comparisons performed for all the neutron star models that we have constructed. These particular plots are drawn for the model #10 of 
the AU EOS, the characteristics of which are presented in Table II of the supplement of Pappas & Apostolatos (2012). Different curves 
correspond to different metrics, i.e., the two-soliton (solid curve), the Manko et al. solution with the negative root for b (dashed-dotted 
curve) and the one with the positive root (dashed-double-dotted), and finally the Hartle-Thorne (dotted) metric. In most cases shown 
here, the Manko et al. curves are on top of the two-soliton curve. The case of the g-t^ component, where the two-soliton curve is almost 
an order of magnitude below the Manko et al. one for a large interval of radii, is a notable exception. The Hartle- Thome's characteristic 
failure to describe gt^ is also evident. 



ric component right at the pole where for some models it's 
fractional difference is a bit smaller than 1/100). In compar- 
ison to the other two analytic metrics discussed above, we 
see that for the models that a Manko et al. metric can be 
found, that is for the rapidly rotating neutron stars, both 
this metric and the two-soliton metric perform very well 
(actually the Manko et al. solution performs better than 
iBerti fc St crgioulas (2 0o3l had in itially found as was shown 
by IPappas fc Apostolatod l|2012h ) and there are only tiny 
differences between the two-soliton and the Manko et al. 
analytic metrics. The greatest difference between the two- 
soliton and the Manko et al. appears to be in the gt^ compo- 
nent of the metric, where the two-soliton is usually clearly 
better. This was anticipated because the gt(f, component of 
the metric is, as we have shown in section 3, more sensitive 
to the value of Ss, which can be suitably adjusted in the 
two-soliton metric, but not in the Manko et al. solution. For 
the rapidly rotating models the Hartle-Thorne is not such 
a good representation of the numerical metric as the other 
two metrics. That also was expected since the Hartle-Thorne 
metric is not suitable for fast rotation. Hartle- Thome's fail- 
ure is more evident in the gt(f, component of the metric which 
is consistent with Hartle-Thorne's vanishing spin octupole 
S-i, (in the Appendix we show that Sz = Q for the Hartle- 
Thorne metric). 

For the slowly rotating models, there are no Manko 
et al. solutions to compare to the numerical metric, so 
in these cases the only alternative is the Hartle-Thorne 
metric. We should say again that the consistent way to 
calculate t he Hartle-Thorne parameters is the one de- 
scribed by iBerti et al.l (|2004l ). but as it is discussed by 



IPappas fc Apostolatod (|2012l ) the parameters of the Hartle- 
Thorne metric (specifically the parameter q which is the re- 
duced quadrupole) consistently calculated are in very good 
agrement with the numerical multipole moments, so these 
moments were used straightforwardly for the construction of 
the Hartle-Thorne metric. Again we saw that the two-soliton 
performs better when compared to the numerical space-time 
than the Hartle-Thorne metric. We should note though that 
the problem of Hartle-Thorne's metric to acurately describe 
the gt4, metric component is present even at slow rotation. 

Having demonstated the overall superiority of two- 
soliton, compared to Manko et al. and Hartle-Thorne, to 
accurately describe the metric functions of any numerical 
neutron star model, in the following comparisons we will 
only compare the two-soliton quantities to the correspond- 
ing numerical ones. 

The next quantity we have used for comparison is the 
position of the ISCO. In Figure [5] we present the relative 
difference between the numerical and the analytic ISCO for 
all neutron star models constructed with the AU EOS for 
both prograde and retrograde orbits (the latter are indicated 
by negative parameter j = J/M"^). The general conclusion 
is that for all models constructed using all three EOSs the 
ISCO of the analytic metric does not deviate by more than 
4 per cent from the ISCO of the corresponding numerical 
model and such deviations are observed for the prograde or- 
bits of the fastest rotating models. We should note that we 
don't perform any comparison with numerical models the 
radius of the surface of which on the equatorial plane is fur- 
ther than the radius of the corresponding ISCO; therefore 
the points that would correspond to these models are miss- 
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Figure 5. These figures show the relative difference between the numerical and the analytic Risco for the neutron star models calculated 
using the AU EOS. The relative difference is shown only for the models that have their ISCO outside the surface of the neutron star. 
The negative values of j correspond to counter-rotating orbits (retrograde). The data for the co-rotating orbits can be found in Table 2. 



ing from the plots. At this point we should mention that 
apart from the position of the marginally stable circular or- 
bit for the particles, there is also the position of the unstable 
photon circular orbit that could also be used as a criterion 
for comparison. This orbit though is usually, for the pro- 
grade case, below the surface of the neutron star (while for 
the retrograde it is usually outside the star) and it doesn't 
have an immediately measurable effeclQ. 

We continue to the results of the comparison of the vari- 
ous frequencies associated to the circular orbits on the equa- 
torial plane, i.e., Q, ^Ip, Hp,^lz, and k^. The analytic orbital 
frequency compares very well to the numerical orbital fre- 
quency for all the models with the relative difference being 
in all cases smaller than ~ 10~^. This result is very impor- 
tant, since the orbital frequency together with the i?isco 
are relevant to observations from accretions discs. Typical 
plots of Q and the relative difference between the analytic 
and the numerical one as functions of the logarithm of the 
distance from the central object are shown in Figure [6] 

With regard to the comparison of the radial and verti- 
cal perturbation frequencies Kp, and Kz, respectively, and the 
corresponding precession frequencies Q,p, and il^, things are 
more complicated. These frequencies include second deriva- 
tives of the metric functions in their calculation. Conse- 
quently the results of the corresponding numerical calcu- 
lations are plagued by accuracy problems. In particular nu- 
merical calculation of the second derivatives induces artifi- 
cial oscillations i n the result s. The se issues have been dis- 
cussed earlier by I Bert i et al.l (|2004l '). Following the sugges- 
tions of N. Stergioulas, we found two ways to mitigate the 
problem. The first one is to calculate the frequencies directly 
in the coordinates that the RNS code produces the metric 
functions so as to avoid any numerical errors caused by the 
transformations of the coordinates and the metric functions 
themselves. Then one would only have to identify the co- 
ordinates of the points at which the frequencies are calcu- 
lated with the corresponding Weyl-Papapetrou coordinates 
(the calculated frequencies themselves are the frequencies 
that static observers at infinity measure so they don't de- 
pend on the coordinate system used). The second one is 
to smooth out these artificial oscillations by taking a three 
point average of the frequencies. The efficiency of this tech- 

^ One could argue that it could be associated to the optics 
around neutron stars and possibly to quasi-normal modes of the 
space-time around the neutron star. 



nique has been tested in the case of the non-rotating mod- 
els (the exterior of which are described by a Schwarzschild 
metric) and has been verified to give trustworthy results. 
Another thing that we should also consider is that in the 
case of Kz, the values are very close to the corresponding 
values of Q; consequently there is low accuracy in the cal- 
culation of the precession frequency fi^ = r2 — in some 
cases. That is why we consider as a better indicator of the 
actual ability of the analytic metric to capture the behavior 
of the numerical metric, the deviations in the oscillation fre- 
quencies Kp and Hz, instead of Qp and ilz- Nevertheless, we 
present both the precession and the oscillation frequencies 
of the analytic metric in comparison to those of the numeri- 
cal metric, together with the relative difference between the 
analytic and the numerical oscillation frequencies, in order 
to get a clearer picture of the comparison. All these plots, 
again for the model #10 of the AU EOS, the characteris- 
tics of which are presented in Table II of the supplement of 
iPappas fc Apostolatod l|2012l '). are shown in Figure [S] 

Generally, the relative differences in Hp between the nu- 
merical and the analytic metrics is small and in some cases 
they could climb up to 10 per cent near the ISCO. This is 
due to the fact that the radial oscillation frequency tends 
to zero as the ISCO is approached, causing an increase of 
the relative difference. In contrast, for Kz, the relative dif- 
ferences are always below 1 per cent at the ISCO. Now the 
picture is inverted for the precession frequencies, which is 
related to the fact that in this case the small quantities are 
the Qz's. The overall picture we obtain is that the analytic 
frequencies capture quite well the behavior of the numeri- 
cal frequencies both qualitatively and quantitatively. This is 
especially evident in the case of the Kz frequency where in 
some cases it becomes greater than the orbital frequency, an 
effect more prominent in the models of EOS L, which the 
two-soliton metric captures quite well. The importance of 
this effect and its relevance to QPOs is further discussed by 
iPappasI l|2012l ). 

The final comparison criterion is the quantity AE. The 
numerical computation of this quantity has similar difficul- 
ties with the precession frequencies; these issues could be 
fixed by performing the same tricks to avoid numerical os- 
cillations. In Figure[7]we show for the same model of a rotat- 
ing neutron star as in the previous cases, the quantity AE 
computed from the numerical and the analytic metric on the 
left, and their relative difference on the right. Again we see 
that the two-soliton metric describes with high accuracy the 
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Figure 6. The plots show the various frequencies and the relative difference of the analytic to the corresponding numerical ones. The 
top row of plots shows on the left the numerical (dotted) and the analytic (solid) orbital frequency Q for comparison reasons (they 
are hardly distinguishable), while on the right it is plotted the logarithm of their relative difference. The middle row shows on the left 
the numerical (dotted) and analytic (solid) precession frequencies Qp, in the middle the corresponding radial oscillations frequencies Kp 
(same correspondence of lines) and on the right the logarithm of the relative difference of the latter oscillation frequencies. Finally the 
bottom row shows the corresponding plots for the plane precession (left) and the vertical oscillation frequencies (middle), Q^, and 
respectively, as well as the logarithm of the relative difference of the latter frequencies (right). All plots are for the model #10 of the AU 
EOS presented by Pappas & Apostolatos (2012). 



AE which we obtain numerically from the numerical mod- 
els. We remind here that this quantity is relevant for the 
emitted spectrum of a thin accretion disc and its tempera- 
ture profile, as well as for the efficiency of the disc, i.e., the 
amount of kinetic energy transformed to radiation. 

We close this section with Table [21 where we present for 
all numerical models of EOS AU (t he multipole moments 
of wh ich are given in the Tables of IPappas fc Apostolatoa 
l|2012l) ) the parameters and the type of the two-soliton met- 
ric along with a few quantities of astrophysical interest, and 
specifically their comparison between the ones calculated us- 
ing the analytic and the numerical space-times. These quan- 
tities are, the circumferential radius at the ISCO i?isco, the 
efficiency rj — 1 — i5isco of a thin accretion disc (if there 
was one around the particular neutron star), the orbital fre- 
quency at the position of the ISCO f2isco (this is a frequency 
expected to show up in QPOs if the latter are related to the 
orbital motion) , and finally the vertical oscillation frequency 
at the ISCO (k2)isco (this could also be related to QPOs). 
The Table shows that the relative differences between the 
numerical and the analytic quantities is of the order of 1 per 
cent or lower. 



6 CONCLUSIONS 

In this work we have tested whether the four-parameter 
two-soliton analytic metric , which was derived by 
iManko. Martin fc Rui j (|l995l ). can be used as a trustwor- 
thy approximation for the space-time around all kind of 
neutron stars. To match the particular analytic metric to 
a specific neutron star mod el, which was produced throug h 
the numerical code RNS of IStergioulas &: FriedmanI l|l995h . 
we have used as matching conditions the first four non-zero 
multipole moments. Our choice was justified, apart from the- 
oretical reasoning, by comparing the numerical metrics with 
different analytic metrics produced by slightly varying two 
of their moments (quadrupole and spin-octupole) . The com- 
parison showed clearly that the best matching comes from 
imposing the condition that the parameters of the analytic 
metric should be such that the analytic space-time acquires 
the first four non-vanishing moments of the numerical met- 
ric. 

Having demonstrated the appropriateness of the match- 
ing conditions, we proceeded to compare the various numeri- 
cal neutron star space-times with the corresponding analytic 
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Table 2. For each neutron star model, that was constructed by the RNS code using the AU EOS, wc have computed a number 
of parameters that are related to the particular analytic two-soliton metric which better approximates the numerical metric. These 
parameters are (i) the type (Case) of the two-soliton space-time that each particular model corresponds to, (ii) the parameters k and b of 
the two-soliton as well as the spin parameter j and, (iii) the relative difference between the following analytic and numerical quantities: 
■RiSCOi JJi ^^ISCO and (kz)iscOi for the models that the ISCO lies outside the surface of the star. The models shown here are the same 
models presented in Pappas & Apostolatos (2012). The rest of the physical parameters of these models, such as the mass and all the 
other multipole moments of the rnotlcls. can be found there. All the lolativc differonces are sivon as a y)erceutago. 



model 


Case 


j 


b 


k 


AKisco 


Ar; 


Ariisco 


(Ak2)isco 








(km) 


{krn'^) 


(per cent) 


(per cent) 


(per cent) 


(per cent) 


1 


III 


0. 


-3. 





0.003 


0.032 


0.006 




2 


la 


0.2015 


-0.0784 


-0.5271 


0.013 


0.019 


0.003 


0.024 


3 


Ila 


0.3126 


-0.1305 


-1.2503 


0.283 


0.142 


0.027 


0.141 


4 


Ila 


0.414 


-0.1858 


-2.1664 










5 


Ila 


0.4749 


-0.224 


-2.8362 










6 


Ila 


0.5297 


-0.2626 


-3.517 










7 


Ila 


0.5789 


-0.3014 


-4.1972 










8 


Ila 


0.617 


-0.3351 


-4.7738 










9 


Ila 


0.651 


-0.368 


-5.3277 










10 


Ila 


0.6618 


-0.3793 


-5.5114 










11 


III 


0. 


-3. 





0.003 


-0.021 


-0.004 




12 


III 


0.194524 


-0.0915 


-0.1528 


-0.021 


0.022 


0.002 


-0.034 


13 


III 


0.309849 


-0.1669 


-0.4288 


-0.018 


0.025 


0.003 


-0.071 


14 


III 


0.406932 


-0.2417 


-0.8163 


0.019 


0.03 


0.005 


0.113 


15 


III 


0.485572 


-0.3158 


-1.2698 


0.148 


0.114 


0.019 


0.04 


16 


III 


0.550214 


-0.3897 


-1.7658 


0.391 


0.204 


0.035 


0.217 


17 


III 


0.603381 


-0.4628 


-2.279 


0.778 


0.463 


0.077 


0.477 


18 


III 


0.645447 


-0.5317 


-2.7705 


1.274 


0.83 


0.137 


0.815 


19 


III 


0.676639 


-0.5916 


-3.197 


1.798 


1.205 


0.197 


1.15 


20 


III 


0.706299 


-0.6585 


-3.6602 


2.453 


1.831 


0.297 


1.638 


21 


III 


0.510282 


-0.3726 


-0.7958 


0.043 


0.026 


0.005 


0.15 


22 


III 


0.510617 


-0.3717 


-0.8179 


0.054 


0.046 


0.008 


0.058 


23 


III 


0.514032 


-0.3731 


-0.8721 


0.063 


0.05 


0.009 


0.031 


24 


III 


0.520506 


-0.3789 


-0.9409 


0.083 


0.016 


0.004 


0.065 


25 


III 


0.547452 


-0.4083 


-1.1827 


0.164 


0.072 


0.013 


0.153 


26 


III 


0.587439 


-0.4607 


-1.5504 


0.346 


0.206 


0.034 


0.164 


27 


III 


0.626593 


-0.5214 


-1.9574 


0.63 


0.391 


0.064 


0.362 


28 


III 


0.659098 


-0.5806 


-2.3502 


0.991 


0.645 


0.104 


0.64 


29 


III 


0.694585 


-0.6577 


-2.8456 


1.551 


1.018 


0.162 


0.948 


30 


III 


0.713165 


-0.7054 


-3.1406 


1.95 


1.349 


0.214 


1.221 
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Figure 8. Plot showing the numerical (dotted) and analytic 
(solid) precession frequencies for the model #10 of the L EOS. 
One can see how the frequencies that are calculated from the ana- 
lytic metric capture the qualitative behavior of the corresponding 
numerical frequencies. The parameters (i.e., the multipole mo- 
ments) for this model can be found in Pappas & Apostolatos 
(2012). 



space-times. To perform the comparison, we have assumed 
several criteria having in mind that they should correspond 
to geometric and physical properties of the space-time, with 
a special interest in physical quantities that could be asso- 
ciated to astrophysical processes that are usually observed 
from the vicinity of neutron stars. 

The result of these comparisons is that the two-soliton 
space-time can reproduce the properties of the space-time 
around realistic neutron stars, and in particular it can re- 
produce all astrophysically interesting properties. Probably 
the most important fact is that the analytic metric can cap- 
ture properties of the neutron star space-time that a corre- 
sponding Kerr space-time could not, such as the behavior 
of the precession frequencies of almost circular and almost 
equatorial orbits. A typical example is shown in Figure 8, 
where we present the analytic and numerical frequencies of 
the precession of the orbital plane for a model constructed 
using the L EOS. The possible importance and implication 
of this, i.e,. the capability of the two-soliton to capture this 
particular behavior in contrast to the Kerr geometry, was 
further discussed in iPappasI (|2012l ). 

Generally the two-soliton metric can be a very useful 
tool for studying phenomena that happen around all kind 
of neutron stars and are quite sensitive to more realistic and 
accurate geometries than the ones used so far. Relying on a 
single analytic metric for all neutron stars is practically more 
favorable than using numerical space-times, or more than 
one analytic metrics depending on the type of the neutron 
star. Thus, the two-soliton metric can be further used for 
more elaborate applic ations such as the ones explored b y 
Psaltis (2008); Baubo ck Psaltis. Ozel fc Johannseiil l|2012l) : 
Baubock. Psaltis fc Ozell l|2012l) . 
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APPENDIX A: NEUTRON STAR MODELS 

In order to construct the analytic space-time exterior to a 
compact object, one has to choose the appropriate mul- 
tipole moments. For neutron stars, these moments could 
be computed from numerical models that are constructed 
with realistic equations of state. There are several schemes 
developed for numerically integrating stellar models (see 
[Stergioulas & Friedman (1995), and for an extended list of 
numerical schemes see IStergioulad (|200^)). We have used 
Stergioulas's RNS code for the construction of the models. 

In order to cover more space on the "neutron star pa- 
rameter space" , we constructed numerical neutron star mod- 
els with equations of state of varying stiffness. For that pur- 
pose we have chosen AU as a typical soft EOS, FPS as a rep- 
resentative moderate stiff EOS, and L to describe stiff EOS. 
Having the numerical models ready, we proceeded in eval- 
uating their multipole mome nts according to the algorithm 
described in lPappas fc Apostolatos ( 2012 ). The parameters 
then used to construct the analytic space-time models, i.e., 
M, a, k, and b, are evaluated from the first four multipole 
moments {M, J, M2, S3) of each model by inverting equa- 
tions (fTU)) . 

For the spe cifics of the various models chosen here, 
we have followed iBerti fc Stergioulas! (|2004l ) . We have con- 
structed the same constant rest-mass sequences as the ones 
presented in lBerti fc StergFo ulas (2004) for the correspond- 
ing equations of state. For every equation of state, 3 se- 
quences of 10 models were constructed, that corresponded 
to: 

(i) a sequence corresponding to a neutron star of 1.4Mq 
in the non-rotating limit, 

(ii) a sequence terminating at the maximum-mass neu- 
tron star in the non-rotating limit, 

(iii) a supermassive sequence that does not terminate at 
a non-rotating model at its lower rotation limit. 

All the sequences end at the mass-shedding limit on 
the side of fast rotation, i.e., at the limit were the angular 
velocity of a particle at the equator is equal to the Keplerian 
velocity at that radius. These sequences are the so called 
evolutionary sequences. 

All the parameters for the comp uted models can be 
found in lPappas fc Apostolatod (|2012h . 



APPENDIX B: THE FUNCTIONS OF THE 
TWO-SOLITON 

In this Appendix we present the full expressions for writing 
the metric functions of the two-soliton. The determinants 



that are given in Sec. 2 and appear in the formulas for the 
metric functions, can be substituted with the following ex- 
pressions, starting with E± — B, where the functions 
A, B are given as: 

A = -16dfc (r_r+ + R-R+) 

- [Ai {R-r+ + r-R+) - (i?_r+ - r_i?+)] 4 

+ [A+ {r.R. + r+R+) - A+ {r.R. - r+R+)] kI, (B1) 

At = (4(a - bf{ab ±d-k)- 4((a - 2b)b ± d)M^) , (B2) 
A2 = ■ii{{a~b){ab±d- k) -bM^) n±, (B3) 

B = 2K-K+M [{R+ - R-)B^ + (r_ - r+) Bf 

+ + R+) B2 - {r- + r+) B+] , (B4) 

Bt = i {2k{a -b) + b (M^ - a'^ + b"^ ± k-k+)) 

X (k+±k_), (B5) 
B2 ^ 2d{M^ -a^ + b'^ ±K-K+) . (B6) 

The determinant H can be substituted as H — —L, where 
L is given by the expression, 

L = {r-r+ + R-R+) Li + (r_ + R- + r+ + R+) L2 

- {r- -R-+r+- R+) L3 + {r-R- - R+r+) 
+ [r- -r+ + R-- R+)Ll + (r__R+ - R-r+) 
+ (r_ -r+-R-+ R+)Li + {R-r+ + r-R+) 

+ {R^r^+r+R+)Le, (B7) 

Li = -16dfc(ia - ib + M - z)M'^, (B8) 
La = i4dM{a + iM) [M* -~2{a -b^ - 2k) 

+ {a-bf {{a + bf-4k)] , (B9) 
L3 = 4dM [ia^ - Ma - i {b^ + + 2k) a + 

+ 2ibk + b'^M] (BIO) 

= 4K=pM {ifeM^ + abM^ 

- i {2ba^ -~ka~b{b^ ±d + 2k)) 

- {2ba^ - ka'^ -b{b^ ±d + 3k) a 
+ fc(2fe^=Fd + fc)) 

+ (2fc^ - (a^ + 3ba T 2d) k + b{a + b) (a^ =F d)) 

X i{a - b)M 

+ {a-b){ab-k) ({a + b) (a"^ T d) + (& - 3a)fc) } (Bll) 

Lf = +4i [b{M - z)a'^ +{b^^d-k){M - z)a^ 
+ ({z- M)b^ +iM^b'^ 
+ (2M^ ±d + 5k){z~ M)b =F idM'^) a^ 
+ ((z - M)b'^ (Af^ =F d - 5fc) (z - M)b^ 
~ 4ikM^b +{k±d) (M^ + 4fc) (M - z)) a 

- ib-^M^ - ib'^M^ {M^ Td-2k) 

+ iM'^ (2k{k ±d)± dM^) + b {l? (M^ ±d + k) 

+ {M"^ + {k±d)Af -4k{k±d))) (M-Z)] Kqr,(B12) 
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= ±4 [b(M - z)a-' + {{k±d){z~ M) + ibM^)a'^ 
+ {2{z - M)b^ + 2 {A'f + 2k) {z - M)b 

- ikM^) 

+ {-2iM^b^ + 2{5k±d){M - z)b'^ 

- i {2M^ + 3kM^) b 

+ {{k±2d)M'^ +Ak{k±d)){M - z))a 

+ { (M - - 2 (2fc - Af^) (M - z)b^ 

+ YikM^b'^ + (M'* + 4kM^ - 8k{k ± d)) {M - z)b 

+ ikM^ {M^ + 2{k ± d))} a 

+ ib^M^ - ib^ {3k - 2M^) 

+ ibM^ {M'^ + 3kM'^ - 2k{k ± d)) 

+ (M^ {2k{k ±d)± dKf) + b'^(k ± dj) (z - M) 

+ b^ {4k{k ± d) - (3fc ± 2d)M^) (M - z)] . (B13) 

The determinant G can be expressed as G = —E, where E 
is given by the expressions, 

E = El {r-r+ + R-R+) + E2 (r_ + R- + r+ + R+) 
+ Ez (r_ -R-+r+- R+) + Et{r-R+ - R-r+) 
+ E^ (r--R- - r+R+) + Et {R-r+ + r-R+) 
+ E^ (r--R- + r+R+) + Et {r- - r+ - R- + R+) 
+ Ee{r.~r+ + R--R+), (B14) 

El = 16dk{-ia + ib + M)M'^, (B15) 

E2 = 4dK^K+(ia + M - z)M, (B16) 

£3 = 4d [2jfc(a - 6) + (6^ + - a^) {M - z + ia)] 

X K-K+M, (B17) 

= -4i [foAf"* + i (b^ =F d) Af^ 

- (a + 6) {2ab - b^ ^ d - k) M"^ 

- i (2k'^ + 2 (6^ - 2ab ±d)k 



+ (a 



+ (a - 6) ((a + 6)^ - 4fc) (ab =F - fc)] Kzp AT, (B18) 

Ef = ±4Af [ibA/^ - {ab^d)M^ 

- i (-26^ + 2a b - Skb - ak) 

+ (2fc^ - (a^ + Aba ~ 3fe^ T 2d) ifc 

+ 2(a^ fe^)(afe=F d)) A/^ 

+ i(a - b) ((a - 6)6(a + bf + 2fc^ 

- (a^ + Aba - 36^ =F 2d) fc) Af 

- (a-6)^((a + 6)^-4fc) (a6=Fd-fc)] , (B19) 

E^ = Ak:^M {6a^ - (fc + i6(Af - z))a'' 

- (b^ + (2Af^ ± d + 3fc) & - ifc(A/ - z)) a^ 

+ {i{M - z)b^ + 5kb^ + i {±d + 2 {A'f + k)) 

X {M - z)b + k (M^ ± d + 3k)) a"^ 

+ {{M^ ±d-k)b^ - 3ik{M - z)6^ 

+ (Af* + (3fc ± d)Af^ - 4fc^) b 



- ifc (A/^ + 2{k ± d)) (Af - z)) a 

- (fc =F d)fcAf^ + fe^fc (-2M^ T d + fc) 

- ib^ (M^ ±d){M- z) 

- ib [M^ + {2k ± d)A/^ - 2k{k ± d)) (M ~ z)} , (B20) 
and finaUy the determinant Ko can be expressed as 

Ko = -16dKiKl. (B21) 

In order to get the metric functions one should substitute 
these expressions to the equations (|ll)ll4p . We should note 
here that the above expressions are not equal to the cor- 
responding determinants since a common factor to all the 
determinants, i.e., the quantity efeeL has been simpli- 

fied out of the expressions. 



APPENDIX C: TRANSFORMATION TO 
WEYL-PAPAPETROU COORDINATES FOR A 
GENERAL METRIC 

As we have mentioned, a stationary and axially symmetric 
space-time g^,^{x^) can be cast in the form of the Papa- 
petrou line element ([1} 

ds^ = -f{dt - ujd(t>f + /"^ (e^^(d/9^ -f- dz^) + p^dct)'^) , 

by an appropriate coordinate transformation, where the 
metric functions are functions of the Weyl-Papapetrou co- 
ordinates (p, z). These coordinates as expressed as func- 
tions of the previous coordinates are harmonic conju- 
gate functions. That is, the coordinate p{x^) is defined as 
— {gt<t>)'^ ~gttg^4> and satisfies the Laplace equation while 
z(a;'') is it's harmonic conjugate. Thus the integrability con- 
ditions for z are the Cauchy-Riemann conditions which can 
be used to calculate th at coordinate. In an earlier work (see 
IPappas fc Apostolatod ([2OO8)) we have shown the correct 
way to integrate these conditions when one has to trans- 
form a metric given in quasi-isotropic coordinates 



, 2 2i/ ,,2 , 

ds — — e dt + e 



2V. , 



- iodtf + e^*" (dr^ + r^de^) , (CI) 

which are commonly used in numerical integrations of the 
Einstein field equations, to the Papapetrou form. In that 
case, the quasi-isotropic coordinates, r, 8 enter the metric in 
a way that can be easily cast in a cartesian form and thus 
the Cauchy-Riemann conditions are given from the usual 
expressions 



dz 
dvj 

dz 



dp 



— , 



where we have de fined the cartesian coordinates zu — 
r sin 6 = r y/l — fi^, — r cos 6 = rfi. 

In the case of a metric given in a general form though, 
like the Hartle-Thorne metric, the Cauchy-Riemann condi- 
tions can not be used as given in the previous equations. So, 
one has to calculate the general form of these conditions. 
The general expressions for the Cauchy-Riemann conditions 
can be evaluated from the orthogonality condition that the 
functions p(r, 9) and z{r, 6) must satisfy. Thus, for a metric 
given in the form 
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ds^ — gttdt^ + 2gt^dtd4> + grrdr^ + 



+ g^^d4P, (C2) 



the orthogonality condition between p and z is Vp ■ Vz = 0, 
which gives the expressions for the Cauchy-Riemann condi- 
tions 



dz 
dr 




(C3) 



(C4) 



where we have expressed the conditions in the r, 6 coordi- 
nates and the function p(r, 6) is defined as 



(C5) 



These general Cauchy-Rie mann conditions are used in the 
same way as prescribed in lPappas fc Apostolatoj (|2008h in 
order to evaluate the z coordinate for the Hartle-Thorne 
metric and then compare the latter metric with the numeri- 
cal metric and the two-soliton metric on the axis of symme- 
try. The metric functions given by the metric in the general 
form can be directly associated to the metric functions in the 
Papapetrou line element ([l]) by comparing the correspond- 
ing gtt, gt(p and g^^ values. The only component that needs 
to be evaluated then is the gpp = g^z and it is given by the 
equation 



r-l 27 

gpp = gzz = f e ' 



(C6) 



APPENDIX D: THE HARTLE-THORNE 
METRIC. 



This is a metric produced bv lHartle fc Thornd (|l968l ) as an 
expansion up to order O(e^), where e = Q/Qt is a parameter 
that characterizes the rotation of a star, and corresponds to 
the exte rior spa ce-time of slowly rotating relativistic stars 
(fi, = M/R^ is the Kepler limit) . The components of the 
metric can be found in lBerti et al.l (|2004l l and are given by 
the expressions, 



2A/\ , 2 



2M ^ 
r J 
2M^ 



qF2) , 



gtt = 

Qrr — 

gee = 

g't>4> ^ 

gt^ = 



where 

- 45r^M-Hl5r*^, 



^ 8A/r4(r ~2A/)' 

- 60/ Af^ + 135r^Af - 45r®) 

+ (16A/^ -f- 8rA'/"* - lOr^Af^ - ?,Qr^M + ISr'^ 



(Dl) 

(D2) 

(D3) 
(D4) 

(D5) 



Ao = 



X {r-M), 

15r(r^2Af)(l^3^^)ln(;-|^) 
16A//2 

15(r2-2A^2) (3u2-l)ln(^) 



16Af2 



Fi = -pW + Ai, 

F2 = 5pr^(r - A/) (2A//^ + 6rAf - 3r^) (3u^ - l) - yli 

Gi = p(-72rA/^-3(L-56A/V)u^-f L) -Ai, 

(16A/^ + 8rA4"* - IQr^WP + Ibr^M + ISr^) 
H\ = 



8A//r4 



X {l~?,u)+A2, 



Ho = 



5 (2M^ - 3rM - 3r2) (l - 3^^) 



8Mr 



A. 



with u — cos 6. We will now try to evaluate the first five mul- 
tipole mome nts of the Hartle-Thorne space-time, following 
iRvanI l|l995h . From the metric functions one can calculate 
the rotation frequency of circular equatorial orbits, which 
is given by equation (|32p . Q, can be expressed with respect 
to a new parameter v = (A/il)^''^. This new parameter is 
a function of r, which can in turn be expressed as a func- 
tion of another new parameter, x — (M/r)^^^. Then the 
parameter v can be expanded as a series in x. This series 
can then be inverted and thus the parameter x can be ex- 
pressed as a power series in v. Similarly, we can calculate the 
energy per mass, of a test particle in circular orbit (|37|l 
as a function of the parameter x. This quantity can then be 
expressed as an expansion on x which by substituting the 
previously obtained power series can then be expressed as 
a series expansion in v. From that, one can calculate the 
invariant quantity 



dn 



V dE 
3 dv ' 



(D6) 



which is related to t he multipole moments of the space-time 
as they are given bv lRvanl (Il995l ). After following all the pre- 
viously mentioned calculations, the expansion of the above 
expression is, 

5 



AP 20jV / 27\ 6 



28jv' 



3 2 
f 80f 



9 



\^ 27 ^ 9 



_ „B ^ _ 81, ^ 



16 



-6gj 



2 J 



( 115 f 35q^ 935g 6615 \ lo 
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\ 18 

1408i^ 
243 



12 24 
572qj 



128 J 



27 



165j v'' +0(v''). (D7) 



From comparing Ryan's expressions to this expansion, we 
see that the Hartle-Thorne metric, as given by Berti, has a 
rotation parameter j defined as j = —J/M^, a quadrupole 
parameter q = Q/M^ and the following moments S3, M4 are 
equal to zero. The above expansion seems to be consistent 
with the aforementioned multipole moments up to 0{v^^). 
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